AD-A255  594 

lilllUlli  _ 

TECHNICAL  REPORT  BRL-TR-3391 


AN  INVERSE  ANALYSIS  OF 
ELECTROTHERMAL-CHEMICAL  (ETC)  GUN 
EXPERIMENTAL  DATA 


GLORIA  P.  WREN 
WILLIAM  F.  OBERLE 


SEPTEMBER  1992 


APPROVED  FOR  PUBUC  RELEASE;  DISTRIBUTION  IS  UNLIMTTED. 


92  9  14  0G9 


iiiiiiiiif 

U.S.  ARMY  LABORATORY  COMMAND 


BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MARYLAND 


NOTICES 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical  Information  Service, 
U.S.  Department  of  Commerce,  5285  Port  Royal  Road,  Springfield,  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position, 
unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers’  names  in  this  report  does  not  constitute  indorsement 
of  any  comm’ercial  product. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

Public  rvponing  burden  for  thi)  collection  of  information  t\  estimated  to  awera^  i  nour  oer  response,  including  the  time  for  reviewing  instruaions,  searching  enstiisg  data  sources, 
gatheemg  end  maintamirtg  the  data  needed,  and  compteting  artd  reviewir>g  the  coiiectton  of  information  Send  comments  reoardirtg  this  burden  estimate  or  any  other  aspect  of  this 
colieaion  of  information,  including  suggestions  for  reducing  this  burden,  to  Washirsgton  neadduarters  Services.  Directorate  for  information  Operations  artd  Reports,  121S  Jefferson 
Oavis  Highway.  Suite  1204.  Arlington.  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-01S8).  Washington.  DC  20S03 

1.  AGENCY  USE  ONLY  (Le»\fe  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 

Sept.  1992  Final,  May  1990  -  Jan  1991 

A.  TITLE  AND  SUBTITLE 

An  Inverse  Analysis  (tf^Etectrothennal-Cheniical  (ETC)  Gun  Eiqierimental  Data 

5.  FUNDING  NUMBERS 

DA311380 

PR:  1F2Z9W9XDGS3 

«.  AUTHOR(S) 

Gloria  P.  Wien  and  William  F.  Obeile 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS{ES) 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

U.S.  Anny  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground.  MD  2100S-S066 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

BRL-TR-3391 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

An  inverse  analysis  is  presented  as  one  method  of  assessing  ballistic  control  exhibited  in  electrothermal-chemical 
(ETC)  experimental  gun  firings.  The  analysis  attempts  to  determine  the  decomposition  or  gas  generation  rate 
necessary  to  satisfy  an  energy  balance  equation  as  gleaned  fiom  available  experimental  data.  The  apparent  functional 
relationsh^  between  the  electrical  energy  input  and  working  fluid  mass  consumed  is  examined.  As  applied  in  this 
report,  the  inverse  analysis  of  ETC  gun  firing  data  is  essentially  an  energy  balance  between  projectile  kinetic  energy, 
gas  kinetic  energy,  and  gas  internal  energy.  Lagrange  assumptions  are  made  for  the  pressure  gradient  and  gas 
velocity.  Thomochemistry  is  assumed  to  be  variable  depending  on  the  plasma  energy  to  winking  fluid  mass  ratio. 

A  sample  database  consisting  of  20-mm  and  30-mm  gun  firings  is  examined. 

14.  SUBJECT  TERMS 

interior  ballistics;  plasma;  guns;  electrical  energy 

15.  NUMBER  OF  PAGES 

81 

16.  PRICE  CODE 

17.  security  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRAa 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UL 

NSN  7S4O-01-280  SS00 


Standard  Form  298  (Rev  2-89) 

^r«Knb«d  bv  ANSI  Sid  ZM-lt 

29B-102 


Intentionally  left  blank. 


TABLE  OF  CO^’''W^S 


Page 


LIST  OF  HGURES  .  v 

ACKNOWLEDGMENTS  .  vii 

1.  INTRODUCTION  .  1 

2.  MODEL  DESCRIPTION  .  2 

3.  NUMERICAL  APPROACH  .  6 

4.  VALIDATION .  7 

5.  THERMOCHEMISTRY  .  11 

6.  RESULTS  OF  INVERSE  ANALYSIS  .  14 

7.  SURFACE  AREA  RATIO .  18 

8.  CONCLUSIONS  .  19 

9.  REFERENCES .  23 

APPENDIX  A:  SAMPLE  INPUT .  25 

APPENDIX  B:  INVERSE  SIMULATION  OF  20-MM  AND  30-MM  FIRINGS 

OF  TiHj/HjO .  29 

APPENDIX  C:  SOURCE  CODE  LISITNG  .  41 

DISTRIBUTION  LIST .  63 


iii 


rrrc  qtt/jjtt  hispectbd  3 


Intentionally  left  blank. 


LIST  OF  FIGURES 


Figure  Page 

I .  Schematic  of  ETC  Gun .  1 

2a.  Comparison  of  IBHVG2  Results  With  the  Inverse  Analysis .  8 

2b.  Comparison  of  IBHVG3  Results  With  the  Inverse  Analysis .  9 

3.  Electrical/Plasma  Energy  for  20-mm  Shot  249  .  10 

4.  Recorded  Pressure  Histories  for  20-mm  Shot  249  .  10 

5.  Comparison  Between  Chamber  Pressure  Gage  History  and  Filtered  Chamber 

Pressure  Gage  History  for  20-mm  Shot  249  .  12 

6.  Projectile  Position  Derived  From  a  Polynomial  Least-Squares  Fit  to  Barrel  Gage 

Arrival  Times  for  20-mm  Shot  249  .  12 

7.  Projectile  Position  Derived  From  a  Polynomial  Least-Squares  Fit  to  Barrel  Gage 

Arrival  Times  for  20-mm  Shot  249  After  133  ps,  Together  With  a  Quadratic 

Fit  From  73  ps  to  133  ps .  13 

8.  Impetus  of  TiH2/H20  With  Electrical  Energy  Added .  14 

9  Electrical  Energy  Input  (kJ)  as  a  Function  of  Time  (ms)  for  20-mm 

Gun  Firings .  16 

10.  Electrical  Energy  Input  (kJ)  as  a  Function  of  Time  (ms)  for  30-mm 

Gun  Firings .  16 

I I .  Normalized  Electrical  Energy  Input  vs.  Nonmalized  Working  Fluid  Consumed  for 

20-mm  Gun  Firings  .  17 

12.  Normalized  Electrical  Energy  Input  vs.  Normalized  Working  Fluid  Consumed  for 

30-mm  Gun  Firings  .  17 

1 3.  Surface  Area  Ratio  for  20-mm  Shots  240,  249,  and  259  Compared  to  Inhibited, 

Monolithic  Grain  .  20 


V 


iNTEmiONALLY  LEFT  BLANK. 


vi 


ACKNOWLEDGMENTS 


The  authors  would  like  to  express  their  appreciation  to  Dr.  Bob  Greig  of  GT-Devices,  a  subsidiary 
of  General  Dynamics  Land  Systems,  who  supplied  the  experimental  data  along  with  needed 
documentation.  Dr.  Paul  Gough’s  criticisms  of  this  work  are  appreciated.  Dr.  Eli  Freedman  of 
Eli  Freedman  Associates  and  Dr.  Kevin  White  of  the  U.S.  Army  Ballistic  Research  Laboratory  (BRL), 
Aberdeen  Proving  Ground,  MD,  supported  this  work  with  thermochemical  data  and  advice.  The  authors 
acknowledge  the  invaluable  assistance  of  Ms.  Sharon  Richardson  of  BRL  who  handled  data  and  performed 
the  code  runs.  Finally,  the  authors  wish  to  thank  Mr.  Fred  Robbins  and  Dr.  Terry  Coffee,  BRL,  for 
reviewing  the  final  manuscript. 


1 


Intentionally  left  blank. 


VIII 


1.  INTRODUCTION 


The  electiotheimal-chemical  (ETC)  gun.  generically  shown  in  Figure  1,  is  a  propulsion  cxincept  which 
utilizes  a  low  mass,  high  energy  plasma  to  initiate  and,  hopefully,  control  the  combustionA'aporization  of 
the  working  fluid  (propellant)  during  the  ballistic  cycle.  Controlling  combustion  (exothermic  working 
fluid)  or  vaporization  (endothermic  working  fluid)  of  the  working  fluid  in  the  ETC  gun  is  necessary  in 
order  to  tailor  the  pressure-time  profile  in  the  gun.  Theoretically,  tailoring  the  pressure-time  profile  to 
obtain  a  "flatter"  and  extended  pressure  curve  slrould  result  in  enhanced  performance  (increased  muzzle 
energy)  and  provide  a  "softer"  launch  environment  to  enhance  projectile  integrity,  especially  for  "smart" 
projectiles.  In  addition,  precise  control  of  the  pressure  history  should  also  allow  for  the  required  projectile 
velocity  repeatability  required  of  indirect  fire  support  (artillery)  applications. 


One  method  of  assessing  the  control  exerted  by  the  plasma  on  the  combustion/  vaporization  process 
is  through  the  use  of  an  inverse  ar  olysis  which  computes  the  decomposition  or  gas  generation  rate  of  the 
working  fluid  given  experimental  measurements  of  pressure,  projectile  motion,  and  plasma  energy  histories 
as  an  input  to  the  model.  If  the  plasma  is  controlling  the  process,  then  the  inverse  analysis  should  show 
decomposition  rates  to  be  a  "strong"  function  of  the  plasma  input  (plasma  energy-time  input). 
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To  illustrate  the  method,  an  inverse  analysis  of  20-mm  and  30-mm  ETC  firings  using  an  exothermic 
working  fluid  is  presented.  The  objective  is  to  ascertain  to  what  extent  these  experimental  ETC  firings 
indicate  control  of  the  working  fluid  combustion  by  the  plasma  energy  input  to  the  gua  The  results  of 
the  inverse  analysis  for  these  data  sets  are  viewed  as  qualitative  since  tlie  experimental  measurements  are 
not  precise;  for  example,  the  projectile  position  is  not  accurately  known  from  the  experimental  data. 

2.  MODEL  DESCRIPTION 

The  derivation  is  traditional  (Hunt  1951).  At  any  given  time,  let  (2  be  the  total  energy  available, 
U  the  internal  energy  of  the  gas,  and  W  the  woik  performed  by  the  gas.  Then 

Q  =  U  +  W  +  Losses  .  (1) 

A  number  of  estimates  exist  for  calculating  losses  due  to  heat  loss,  recoil,  projectile  resistance, 
projectile  rotational  energy,  shocked  air  ahead  of  the  projectile  and  energy  in  compression  of  the  working 
fluid.  In  small  caliber  guns,  heat  loss  can  be  significant.  However,  in  the  current  implementation  for  a 
qualitative  comparison  between  a  number  of  small  caliber  shots,  losses  are  assumed  to  be  similar  in  each 
firing.  Thus,  the  "losses"  term  is  taken  to  be  zero. 

The  total  available  energy  at  the  time  in  which  a  fraction  z(t)  of  a  total  charge  mass  C  has  been  burnt 
in  a  closed  vessel  is 

Q  =  Cz(0  r"  c,{T)dT 

where  c^,  is  the  specific  heat  at  constant  volume  of  the  gases  at  any  temperature  T,  is  the  flame 
temperature  and  is  a  reference  temperature.  Since  only  changes  of  energy  will  be  considered,  the 
actual  value  for  is  not  significant. 

Since  work  is  performed  in  the  gun,  the  gas  temperature  T  is  less  than  the  flame  temperature  T^,  and 
the  internal  energy  is  given  by 
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u  «  CzU)  [J  c^(T)dT  . 

ref 


(3) 


The  woik  performed  by  the  gas  is  the  kinetic  energy  of  the  projectile  and  kinetic  energy  of  the  gas. 
The  kinetic  energy  of  the  projectile  at  any  time  is 


KE(projectile) 


1  2 .  . 
-mvpixj) 


(4) 


where  Vp(x,t)  is  the  projectile  velocity  and  m  is  the  projectile  mass. 

The  kinetic  energy  of  the  gas,  KE(gas).  is  determined  in  a  similar  manner.  However,  since  both  the 
velocity  and  mass  of  gas  can  vary  along  the  length  of  the  gun, 

KE(gas)  =  ^j^P^‘'>pix,t)  A(x)v^(x,t)dx  (5) 


where  x^  is  the  projectile  base  position  at  time  t,  p(x,t),  Vg(x,t)  and  A(x)  are  respectively  the  fluid  density, 
fluid  velocity,  and  cross  sectional  area  of  the  gun  at  position  x. 


To  integrate  Equation  5,  the  Lagrange  assumption  that  the  gas  velocity  Vg(x,t)  is  linear  from  the 
breech  at  position  0.0  to  the  projectile  base  at  x^.  with  constant  spatial  density  p(t),  that  is. 


Vg{x,t)  ^  JLvpit)  and 


ap(jr,r) 

■~5r~ 


=  0 


(6) 


is  utilized.  Then,  assuming  a  constant  area  A  for  the  gun.  Equation  5,  upon  substitution  of  Equation  6, 
becomes 


KE(gas)  « 

Jo  2 


— —  v„(r) 


dx 


J 


(7) 
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(8) 


^  Xp(t) 


L 


or. 


=  -ipCO  A  Xp(t)Vp(t) 


(9) 


KEigas)  *  ^Cv^(t) 


(10) 


assuming  that  the  total  charge  is  accelerated  and  uniformly  distributed.  Substituting  Equations  2,  3, 4,  and 
10  into  Equation  1, 

Cr(r)  C^c.mdT  -  CzO)  T  c^(T)dT  *  Imv^r)  +  lcv^(0  (il) 


Cr(r)  c,{T)dT  =  Imv^^r)  *  jCv^(0  •  '^2) 

Assuming  an  average  value  of  in  the  temperature  range,  T(t)  to  yields 

Cz{t)-,{T^-nt))  -  Imv^^r)  +  -icv^^O  .  (13) 


If  Vp(t)  can  be  obtained  experimentally,  then  Equation  13  has  two  rnknowns,  z  and  T.  A  Nobel- Able 
covolume  equation  of  state. 
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(14) 


P{t) 


V  +  X 


it)A-  —  -  _  bCz  U) 

Pa  Pc(0 


RT(t) 

MW 


CzU) 


provides  a  second  condition  on  the  system.  In  Equation  14,  P(r)  represents  the  space-mean  pressure 
based  on  the  Lagrange  ballistic  assumptions.  the  original  free  chamber  volume,  Xp(t)  the  projectile 
travel,  A  the  cross-sectional  area  of  the  tube,  the  mass  of  air  in  the  chamber,  the  density  of  air,  p^(t) 
the  woridng  fluid  or  charge  density  which  has  been  calculated  from  a  compressible  equation  of  state,  b 
the  covolume  correaion,  R  the  universal  gas  constant,  and  MW  the  molecular  weight  of  reaction  products. 
The  effect  of  air  in  the  chamber,  /p^,  is  included  since  ETC  firings  can  have  large  initial  ullage. 


The  equation  of  state  for  the  unvaporized  fluid  is  in  the  form  (Gough  1989) 


Pc(0 


1  + 


(P(t)-Po) 


//  AT-  =  0 


Po 


1  *!^(Pit)-P^) 


if  K2*0 


(15) 


where  P(t)  is  the  breech  pressure  and  p^  is  the  fluid  density  at  the  initial  breech  pressure  P^  (ambient  in 
these  calculations)  assuming  that  the  bulk  modulus,  K^,  is  at  most  a  linear  function  of  pressure. 


K  +  K2  (P(t)  -  P^  . 


(16) 


RT„ 

Using  the  standard  definition  of  impetus,  /  « _ ,  Equation  14  is  solved  for  T. 

MW 


T(t) 


T7FU) 


V^*x  {t)A-!H-  -bCz(t) 

^  Pa  Pc(') 


(17) 


lCz(t) 
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The  average  value  of  is  expressed  as 


/ 


To(y-i) 


(18) 


where  7  is  the  ratio  of  specific  heats. 

Substituting  Equations  17  and  18  into  Equation  13  and  solving  for  z(t)  at  any  time  yields 


2(0  = 


-mvp(r) 

2  P 


f 

\ 

f  1  c  'I 

1  + 1  — 

[  3  mj 

^  Pit) 

tfi„ 

V  +  X  it) A  -_f. 

C 

(7-1) 

1  '  p« 

Pc(0  J 

1 


(7-1) 


/ 

(  \\ 

C/  +  Pit) 

bC-  ^ 

9  At) 

\  ^  J ) 

(19) 


Experimental  data  are  utilized  to  solve  Equation  19  directly  for  2  by  providing  chamber  pressure  from 
which  the  space-mean  pressure  P(t)  is  calculated  assuming  a  Lagrange  gradient.  Experimental  projectile 
position,  Xp(t},  is  utilized  not  only  to  calculate  the  projectile  velocity  Vp((},  but  also  to  detennine  the  gun 
volume  at  any  time.  The  initial  charge  mass,  C,  the  bore  area.  A,  and  the  projectile  mass,  m,  are 
determined  from  the  experiment.  Propellant  density,  pc(t).  is  calculated  from  Equation  15.  It  is  noted  that 
in  the  calculation,  each  timestep  is  independent  of  any  other  timestep.  The  specification  of  the 
thermochemical  values  (impetus  /,  ratio  of  specific  heats  7,  and  covolume  b)  is  discussed  in  detail  below. 

3.  NUMERICAL  APPROACH 

As  mentioned  earlier,  each  time  step  is  independent  of  any  other  time  step.  Thus,  Equation  19  can 
be  solved  algebraically  for  2(1).  However,  since  impetus,  /,  ratio  of  specific  heat,  7,  and  covolume,  b,  are 
functions  of  z(t)  (see  Section  5),  it  is  necessary  to  iterate  several  times.  Basically,  to  compute  r„+/(/)  at 
time  values  for  /,  7,  and  b  at  time  are  used  in  Equation  19  to  determine  New  values  of 

/,  7  and  b  are  then  determined  based  upon  and  used  in  Equation  19  to  determine  the  next 

approximation.  to  z„_^,j(t).  This  process  continues  until  successive  approximations  are  within  the 

derived  tolerance. 
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4.  VALIDATION 


In  order  to  validate  the  inverse  analysis,  a  120-mm  gun  firing  was  simulated  with  a  standard  lumped 
parameter  interior  ballistic  code,  IBHVG2  (Anderson  1987)  developed  by  the  U.S.  Army  Ballistic 
Research  Laboratory  (BRL).  The  resulting  IBHVG2  breech  pressure  history  and  projectile  position  were 
used  as  input  to  the  inverse  code.  A  comparison  of  the  IBHVG2  calculation  and  the  inverse  analysis 
(ETINV)  for  mass  consumed  is  shown  in  Figure  2.  As  can  be  seen,  the  inverse  analysis  accurately 
reproduces  the  time  history  of  propellant  mass  burnt  However,  as  discussed  elsewhere  (Gough  1990), 
inaccuracies  in  the  experimental  data  can  lead  to  variability  in  the  amount  of  mass  burnt  as  determined 
by  an  inverse  analysis.  Thus,  as  discussed  above,  the  inverse  analysis  is  best  viewed  as  a  qualitative 
picture  of  the  combustion  process. 

The  ETC  experimental  database  examined  in  this  report  consists  of  six  20-mm  shots  using  titanium 
hydride  and  water,  and  three  30-mm  shots  using  titanium  hydride  and  water.  All  firings  were  performed 
at  GT- Devices  (Greig  1989),  a  subsidiary  of  General  Dynamics  Land  Systems.  A  summary  of  the  major 
physical  parameters  for  the  shots  analyzed  is  shown  in  Table  1. 


Table  1.  Summary  of  21 -mm  and  30-mm  Initial  Conditions 


Tube 

Diameter 

Chamber 

Volume 

Working  Fluid 

Projectile 

Mass 

Working  Fluid 
Mass 

Projectile 

Travel 

Maximum 
EE  Input 

21  mm 

72.4  cm^ 

50  %  TiHj  &  50%  HjO 

20  g 

64.0  to  64.64  g 

199  cm 

-500  kJ 

30  mm 

206.0  cm^ 

50%  TiHj  &  50%  HjO 

100  g 

168.84  g 

300  cm 

-400  kJ 

A  number  of  experimental  measurements  is  available  for  each  shot.  The  experimental  data  is  used 
to  determine  plasma  energy  delivered  to  the  combustion  chamber,  chamber  pressure,  and  projectile 
position  as  functions  of  time. 

As  an  example,  the  electrical  or  plasma  energy  computed  from  the  product  of  the  experimental 
voltage  and  current  for  a  20-mm  shot.  Shot  249,  is  shown  in  Figure  3.  Pressure  histories  for  Shot  249 
are  shown  in  Figure  4.  Gages  ate  located  in  the  combustion  chamber,  and  in  the  barrel  at  13  cm,  14  cm, 
18  cm,  88  cm,  and  168  cm  of  projectile  travel.  The  initial  position  of  the  projectile  is  at  tube  origin. 
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Mass  Burned  Computed  by  IBHVG2 

120mm.  no  losses 


Time  Ims) 
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Propelani  Burnt  (g) 
(Thousands) 


Mass  Burned  Computed  by  Inverse  Code 

9.0 

8.0 

7.0 

6.0 

5.0 

4.0 

3.0 

2.0 

1.0 

0.0 

0.0  ZO  4.0  6.0  8.0 

Tims  (ms) 


t20mm.  no  losses 


Figure  2b.  Comparison  of  IBHVG2  Results  With  the  Inverse  Analysis. 
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PRESSURE:  PSI  kilo  JOULES 

(Thouaonda) 


For  the  majority  of  firings  analyzed,  the  chamber  pressure  data  is  smoothed  since  the  pressure  history 
is  "noisy."  To  illustrate,  consider  the  chamber  pressure  for  Shot  249,  Figure  4,  which  is  particularly 
"noisy"  for  the  first  47  ps.  To  perform  the  smoothing,  a  low-pass  filter  is  applied  to  the  data,  and  the 
pressure  history  resulting  from  an  inverse  Fast  Fourier  Transform  (FFT)  is  used  in  the  analysis,  a 
technique  commonly  utilized  in  the  liquid  propellant  gun  community.  A  comparison  between  the  original 
data  set  and  the  data  set  resulting  from  the  inverse  FFT  is  shown  in  Figure  5.  The  filtering  process  is  seen 
to  maintain  the  characteristics  of  the  data  while  removing  the  high  frequency  noise. 

Projectile  position  is  derived  from  arrival  times  at  the  barrel  gage  locations.  Using  the  time  of  arrival 
at  the  barrel  gage  location,  together  with  the  actual  gage  locations,  an  interpolating  fourth-degree 
polynomial  of  position  as  a  function  of  time  can  be  obtained.  Experimental  muzzle  velocity,  as  measured 
by  laser-velocity  trap  data,  correlates  very  well  with  velocity  predicted  by  the  least-square  fit  (Greig  1989). 
In  addition,  any  available  radar  data  are  utilized  to  confinn  projectile  position.  However,  in  general,  the 
radar  data  are  not  of  sufficient  fidelity  to  utilize  in  the  inverse  model.  The  polynomial  method  of 
determining  projectile  travel  distorts  the  initial  estimate  of  projectile  motion  prior  to  the  first  pressure 
probe  since  it  is  an  extrapolation,  as  shown  in  Figure  6  for  Shot  249.  The  projectile  begins  motion  at 
73  ps.  and  the  projectile  travel  is  negative  from  73  ps  to  133  ps.  To  eliminate  negative  projectile  travel 
in  the  analysis,  a  quadratic  fit  is  applied  over  the  time  from  73  ps  to  133  ps  using  zero  travel  at  73  ps 
and  the  calculated  travel  at  133  ps.  The  resulting  estimated  projectile  travel  history  is  shown  in  Figure  7 
for  Shot  249.  However,  this  procedure  calls  into  question  results  obtained  for  approximately  the  first  15% 
of  the  analysis  as  a  function  of  time. 

The  accuracy  of  the  inverse  analysis  is  dependent  on  accurate  experimental  data  for  plasma  energy 
delivered  to  the  combustion  chamber,  chamber  pressure  history,  and  projectile  travel.  Thus,  since  much 
of  the  needed  information  from  the  titanium  hydride  and  water  firings  is  inferred,  as  opposed  to  directly 
measured,  the  analysis  of  these  firings  is  useful  only  as  a  comparison  between  experiments  with  similar 
initial  conditions  (i.e.,  gun  geometry,  plasma  energy,  etc.).  An  accurate  estimation  of  gas  generation  rate 
will  require  precise  experimental  data. 

5.  THERMOCHEMISTRY 

The  specification  of  the  thermochemical  properties  of  the  working  fluid  must  consider  the  input  of  the 
plasma.  Therefore,  a  series  of  BLAKE  (Freedman,  private  communication)  runs  are  performed  with 
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Figiire  5.  Comparison  Between  Chamber  Pressure  Gage  History  and  Filtered  Chamber  Pressure  Gage 
History  for  20-mm  Shot  249. 
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Figure  6.  Projectile  Position  Derived  From  a  Polynomial  Least-Squares  Fit  to  Barrel  Gage  Arrival 
Times  for  20-mm  Shot  249. 
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Figure  7.  Projectile  Position  Derived  From  a  Pplvnomial  Least-Squares  Fit  to  Barrel  Gage  Arrival 
Times  for  20-mm  Shot  249  After  133  us.  Together  With  a  Quadratic  Fit  From  73  us  to 
133  us. 

varying  amounts  of  plasma  energy  input  (Bunte  and  Oberie  1989).  In  the  BLAKE  calculations,  the  plasma 
energy  is  assumed  to  mix  uniformly  with  some  amount  of  woildng  fluid,  and  the  ratio  of  electrical  energy 
input  to  working  fluid  mass  is  referred  to  as  the  electrical  energy  density.  Based  on  the  electrical  energy 
density,  BLAKE  then  deflnes  an  impetus,  ratio  of  specifle  heats,  temperature  and  molecular  weight  for 
the  mixture.  In  the  inverse  analysis,  the  plasma  energy  is  assumed  to  mix  only  with  the  dissociated 
working  fluid  as  determined  iteratively  from  Equation  19.  Since  the  BLAKE  entries  are  not  defined  above 
approximately  6,000  K  where  ionization  would  occur,  the  impetus  values  are  extended  with  a  linear 
regression  line  to  define  electrical  energy  densities  outside  the  table  range. 

Impetus  for  a  mixture  of  50%  TiH2  and  50%  H2O  as  a  fimetion  of  electrical  energy  density  is  shown 
in  Figure  8  for  the  BLAKE  calculations  and  is  denoted  by  the  square  boxes.  The  linear  regression  line 
is  shown  in  Figure  8  as  well,  and  is  given  by 

Impetus  =  134238  *  (EE  density)  +  424.69. 
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Figure  8.  Impetus  of  TiH^/H^O  With  Electrical  Energy  Added. 

The  variable  thermochemistry  values  are  shown  in  Table  2.  As  can  be  seen  from  an  examination  of 
Table  2.  the  values  of  gamma  and  the  covolume  for  electrical  energy  densities  producing  gas  temperatures 
above  5,000  K  are  held  constant  at  the  last  BLAKE  calculation  at  9  kJ/g.  Although  not  accurate,  better 
thermochemical  data  does  not  exist  at  the  present  time. 

6.  RESULTS  OF  INVERSE  ANALYSIS 

The  inverse  analysis  assesses  ballistic  control  exhibited  in  the  experiment  by  inferring  the  relationship 
between  electrical  energy  added  to  the  system  by  the  plasma  and  the  amount  of  woildng  fluid  consumed. 
The  delivery  of  electrical  energy  to  the  system  is  fust  discussed  for  the  experimental  shots.  Results  of 
the  inverse  analysis  are  then  shown  for  20-mm  and  30-tnm  gun  data  with  working  fluid  titanium  hydride 
and  water. 

Electrical  energy  input  in  kilojoules  vs.  time  in  miUiseconds  is  shown  in  Figure  9  for  the  six  20-mm 
firings  (Shots  240, 249, 250, 257, 258,  and  259).  The  shape  of  the  curves  is  quite  similar  with  a  variation 
in  pulse  length  and  total  energy.  The  30-mm  experiments  are  firings  in  which  initial  conditions  and  the 
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Table  2.  Thermochemistry  of  50%  TiHj  and  50%  H2O 


EE  Density 
(KJ/g) 


Temperature 

(K) 


1.627 

1.950 

2.162 

2.630 

3.078 

3.502 

3.899 

4.265 

4.604 

4.917 


Impetus 

(J/g) 

Gamma 

(-) 

Covolume 

(cm^/g) 

511 

1.1685 

0.655 

612 

1.1685 

0.655 

679 

1.2722 

0.655 

826 

1.2553 

0.668 

969 

1.2443 

0.675 

1.106 

1.2370 

0.681 

1.239 

1.2325 

0.690 

1.368 

1.2302 

0.702 

1.495 

1.2294 

0.717 

1.622 

1.2299 

0.735 

1,767 

1.2299 

0.735 

3.109 

1.2299 

0.735 

4,451 

1.2299 

0.735 

5.794 

1.2299 

0.735 

7,136 

1.2299 

0.735 

8.478 

1.2299 

0.735 

9.821 

1.2299 

0.735 

11.163 

1.2299 

0.735 

12.506 

1.2299 

0.735 

13,848 

1.2299 

0.735 

delivery  of  electrical  energy  are  nearly  the  same.  Thus,  these  shots  (Shots  39.  40,  and  43A)  form  a 
repeatability  series  for  this  gun  configuration.  As  shown  in  Figure  10,  approximately  350  kJ  of  electrical 
energy  is  delivered  as  a  simple  ramped  input  over  approximately  1.4  ms. 

To  assess  the  dependence  of  working  fluid  consumed  on  plasma  energy  input,  the  results  are  presented 
as  groups  of  normalized  electrical  energy  input  vs.  normalized  working  fluid  consumed.  The  normalized 
curve  (Figures  11  and  12)  is  obtained  by  dividing  electrical  energy  by  the  maximum  electrical  energy 
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Figure  11.  Nonnalized  Electrical  Energy  Input  vs.  Nomialized  Working  Fluid  Consumed  for  20-iiun 
Gun  Firings. 


Figure  12.  Normalized  Electrical  Energy  Input  vs.  Nonnalized  Working  Fluid  Consumed  for  SO-mm 
Gun  Firings. 


input  and  the  mass  consumed  by  the  maximum  mass  of  working  fluid  in  the  experiment.  The  computed 
mass  of  working  fluid  consumed  is  shown  in  Appendix  B  for  each  shot,  with  the  sample  input/output  in 
Appendix  A. 

The  results  for  the  20-mm  shots  with  titanium  hydride  and  water  as  the  working  fluid  are  shown  in 
Figure  11.  From  Figure  9  the  electrical  energy  input  is  a  simple  ramp.  Although  the  mass  of  working 
fluid  consumed  (Figure  11)  does  indicate  some  linearity  as  a  fimction  of  electrical  energy,  the  results 
provide  no  evidence  that  the  plasma  energy  input  drives  the  ballistic  response  for  the  entire  ballistic  cycle. 

The  30-mm  shots  using  titanium  hydride  and  water  as  the  working  fluid  are  shown  in  Figure  12. 
These  shots  form  a  set  of  repeatability  experiments.  As  shown  in  Figure  10.  the  approximately  350  kJ 
of  electrical  energy  is  deposited  uniformly  over  approximately  1.2  ms  in  each  experiment.  As  can  be  seen 
in  Figure  12.  the  working  fluid  response  is  similar  in  all  three  experiments.  As  in  the  20-mm  experiments, 
although  the  plasma  energy  is  applied  uniformly,  the  rate  of  working  fluid  mass  consumed  appears 
toincrease  during  the  ballistic  cycle.  Although  the  mass  fraction  of  working  fluid  consumed  appears  to 
occur  similarly  in  the  experiments,  it  is  not  evident  that  the  plasma  controls  the  ballistic  process. 

In  both  experimental  data  sets,  an  increase  in  the  rate  of  working  fluid  consumed  over  time  is 
observed.  In  order  to  identify  the  role  of  plasma  energy  as  the  driving  mechanism  in  the  gas  generation 
rate,  shots  in  which  the  initial  plasma  input  is  substantially  varied  should  be  flred  and  compared. 

An  analysis  of  several  firings  with  water  as  the  working  fluid  was  attempted.  It  is  possible  to  devise 
an  extended  thermochemical  description  of  the  mixture  of  water  and  electrical  energy  similar  to  Table  2. 
However,  as  demonstrated  by  Gough  (1990)  and  verified  in  the  analysis  of  experimental  data,  the 
dependence  of  pressure  on  the  fraction  of  working  fluid  consumed  is  weak  for  endothermic  fluids  for  all 
electrical  energy  densities.  As  opposed  to  exothermic  fluids,  the  possibility  of  multiple  roots  and 
bifurcation  in  the  inverse  analysis  makes  water  firings  unsuitable  for  inverse  analysis.  Thus,  these  firings 
are  not  further  discussed  in  this  report. 

7.  SURFACE  AREA  RATIO 

A  second  method  to  investigate  control  of  the  ballistic  process  is  to  ascertain  the  surface  area  of  the 
reacting  working  fluid.  If  the  process  is  not  controlled,  it  would  be  expected  that  the  total  reacting  surface 
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area  would  display  large  fluctuations.  To  determine  an  apparent  surface  area,  the  experimental  breech 
pressure  and  working  fluid  mass  consumed  (as  computed  from  the  inverse  analysis)  are  used  to  compute 
a  surface  area  ratio  (current  surface  area  to  initial  surface  area)  under  the  assumption  of  a  pressure 
dependent  combustion  law  for  the  working  fluid  reaction.  Although  the  results  are  qualitative,  the  ratio 
provides  some  insight  into  the  behavior  of  the  surface  which  is  in  contact  with  the  hot  gases.  Surface  area 
ratios  in  solid  propellants  are  well  behaved. 

A  pressure  dependent  bum  rate  of  JA2  propellant  (r  =  0.12  r  in  cm/s)  and  an  initial  surface 
area  equal  to  the  cross-sectional  area  of  the  chamber  are  assumed.  Then  for  each  percent  mass  consumed, 
the  ratio  of  the  current  to  the  initial  surface  area  can  be  computed.  The  selection  of  the  initial  surface  area 
as  the  cross-sectional  area  of  the  chamber  is  arbitrary.  However,  changing  the  value  of  the  initial  surface 
area  preserves  the  shape  of  the  surface  area  ratio  curve  but  changes  its  magnitude.  Thus,  the  surface  area 
ratio  curves  should  not  be  viewed  in  terms  of  absolute  magnitude  but  only  in  terms  of  characteristic  shape. 
The  results  are  similarly  affected  by  the  choice  of  bum  rate  law.  Figure  13  shows  the  surface  area  ratio 
for  three  20-mm  firings.  Shots  240,  249,  and  259,  with  titanium  hydride  and  water  as  the  working  fluid. 
For  the  first  50%  of  the  mass  consumed  the  ratios  show  minor  fluctuations  except  for  the  early  portions 
of  the  graph  where  the  data  is  questionable  due  to  the  inaccuracies  in  the  early  portion  of  projectile  travel. 
Unfortunately,  the  surface  area  ratios  fluctuate  sharply  after  approximately  50%  of  the  working  fluid  mass 
has  been  consumed,  especially  Shot  249.  In  the  experiments,  the  plasma  continues  to  enter  the  chamber 
until  almost  all  of  the  working  fluid  has  been  consumed  as  seen  in  Figure  11.  Thus,  the  working  fluid 
surface  area  appears  to  be  well  behaved  only  initially. 

For  comparison,  the  geometry  of  an  inhibited,  monolithic  grain  with  dimensions  suggested  by  the 
initial  working  fluid  configuration  is  theoretically  burned  under  the  same  pressure-dependent  conditions. 
The  surface  area  ratio  as  a  function  of  percent  mass  consumed  is  shown  in  Figure  13  and  is  seen  to  be 
progressive  and  well-behaved. 

8.  CONCLUSIONS 


An  inverse  analy.sis  has  been  presented  as  one  method  of  assessing  ballistic  control  exhibited  in  ETC 
experimental  gun  firings.  The  analysis  attempts  to  determine  the  decomposition  or  gas  generation  rate 
necessary  to  satisfy  an  energy  balance  equation  as  gleaned  from  available  experimental  data.  The 
quantitative  values  of  woricing  fluid  mass  consumed  are  dependent  on  accurate  experimental  breech 
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Grain. 

pressure  measurements  and  projectile  position  as  well  as  the  assumption  of  the  form  of  the  pressure 
gradient 

The  analysis  was  applied  to  20-mm  and  30-mm  experimental  firings  (rear  plasma  injection)  using 
titanium  hydride  and  water  as  the  working  fluids.  The  results  of  the  analysis  are  viewed  as  qualitative 
for  these  data  sets  and  suggest: 

(1)  For  the  firings  analyzed,  there  is  little  evidence  to  indicate  that  the  electrical  energy  input  plays 
a  dominant  role  in  the  amount  of  working  fluid  consumed. 

(2)  A  large  fluctuation  in  surface  area  occurs  for  the  titanium  hydride  and  water  system  when 
approximately  50%  of  the  mass  has  been  consumed  regardless  of  the  presence  or  absence  of 
plasma  energy  in  the  configuration  analyzed. 
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(3)  Accurate  measurement  of  projectile  position  through  in-bore  radar  is  needed  to  accurately  establish 
the  volume  behind  the  projectile. 

(4)  The  potential  coupling  between  plasma  energy  input  and  working  fluid  mass  consumed  is 
obscured  by  the  lack  of  comparable  data  sets  in  which  the  plasma  energy  input  history  is 
significantly  varied. 

(5)  Additional  work  is  needed  in  defining  the  thermochemistry  of  working  fluid/electrical  energy 
combinations. 
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APPENDIX:  A 
SAMPLE  INPUT 
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Enter  the  file  containing  Input  file  names. 


Inverse  Electrothermal  Gun  Code 


Pressure -time  file:2S0tvp.IN 
ProJ  motion* time  f ile : 250tvpp . IN 
Proj  travel -time  file. 

Electrical  energy- time  file:250tvce.IN 

Initial  data  flle:250inlt.IN 

Chamber  volume  (cm*3)-  72.40000 

Tube  diameter  (cm)-  2.1  Area  (cm*2)-  3.46361 

Projectile  position  (cm)-  0.00000 

Projectile  mass  (kg)-  0.02 

Working  fluid  mass  (g)-  64.64000 

Initial  working  fluid  density  (g/cm*3)-  1.59000 

Tolerance  required  in  mass  consumed-  0.00100 

Time  start  (sec )-0. 00000000 

Time  Increment  (sec)-O. 00005000 

Time  end  (sec)-O . 00120000 

Initial  guess  for  mass  consumed  (g)-  0.10000 

Initial  value  of  beta,  fraction  of  wf  moving  (-)-  1.00000 

Mass  of  air  initially  in  chimber  (g)-  0.03750 

Bulk  modulus  of  WF  at  r  atmosphere  (MPa) -5000. 00000 

Derivative  of  bulk  mod'  1 us  wrt  pressure  (-)-  8.00000 


Thermochemical  dnta  file:tlhextra.prn 


TIME 

CHAM  PR 

BASE  PR 

PROJ  TRAV 

PROJ  VEL 

EL  EN  DEN 

(ms) 

(MPa) 

(MPa) 

(cm) 

(cm/s) 

(MJAg) 

0.00 

0.10 

0.04005 

0.0000 

0.00 

0.000 

0.05 

2.89 

1.15726 

0.3916 

7831.56 

19259.860 

0.10 

98.12 

39.29342 

0.7832 

7831.56 

2101.635 

0.15 

200.20 

80.17240 

1.1747 

7831.56 

1811.112 

0.20 

338.07 

135.38190 

1.5663 

27364.79 

1433.635 

0.25 

413.63 

165.64180 

3.9112 

55843.03 

1198.416 

0.30 

493.50 

197.62530 

7.1506 

73255.77 

1039.150 

TIME- 

0.00035 

NO  CONVERGENCE 

0.40 

573.84 

229.79470 

16.1227 

105250.30 

17.629 

0.45 

577.92 

231.42860 

21.7618 

119858.60 

11.768 

0.50 

588.95 

235.84870 

28.1085 

133558.50 

8.616 

0.55 

559.01 

223.85940 

35.1177 

146363.10 

7.105 

0.60 

473.38 

189.56690 

42.7448 

158286.10 

7.172 

0.65 

403.04 

161.39790 

50.9463 

169340.30 

7.109 

0.70 

355.40 

142.31970 

59.6788 

179538.80 

6.588 

0.75 

306.18 

122.61100 

68.9001 

188895.20 

6.544 

0.80 

263.75 

105.61900 

78.5683 

197422.40 

6.623 

0.85 

235.21 

94.18995 

88.6424 

205134.10 

6.418 

All  fluid  consumed. 

Solution  continues  to  balance  energy. 

0.90 

213.39 

85.45189 

99.0817 

212043.00 

6.418 

0.95 

190.69 

76.36062 

109.8466 

218163.30 

6.418 

1.00 

171.73 

68.76966 

120.8980 

223505.80 

6.417 

1.05 

156.52 

62.67997 

132.1973 

228086.30 

6.417 

1.10 

143.31 

57.38857 

143.7066 

231917.30 

6.417 

1.15 

130.72 

52.34921 

155.3890 

235012.50 

6.416 

1.20 

119.08 

47.68684 

167.2078 

237381.90 

6.416 

CUM  WF 
(g) 

0.00000 

0.00079 

0.02080 

0.04470 

0.08227 

0.12971 

0.19483 

16.52764 

28.05256 

43.34926 

56.13565 

56.85619 

57.92583 

62.86549 

63.38555 

62.64100 

64.93068 


66.28542 

66.49628 

66.69440 

67.08662 

67.33691 

67.08495 

66.46645 
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APPENDIX  B: 

INVERSE  SIMULATION  OF  20-MM  AND  30-MM  FIRINGS  OF  TfHj/HjO 
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APPENDIX  C: 
SOURCE  CODE  LISTING 


41 


Intentionally  left  blank. 


42 


PROGRAM  ETINVBC 
C 

C  ADD  COMPRESSION  OF  WORKING  FLUID  -  JULY  1990 
C 

C  CONSIDER  BETA~AMT  OF  FLUID  WHICH  MOVES 
C  APRIL  1990 

C  INVERSE  PROGRAM  TO  CALCULATE  INTERACnON  RATE 
C  OF  WORKING  FLUID  AND  PLASMA  FROM  ET  DATA 
C  ASSUMING  LAGRANGE  PRESSURE  GRADIENT 
C 

C  TYPE  ETINVBC  (WAIT  FOR  QUESTION  ABOUT  INPUT  FILE!!) 

C 

IMPLICIT  REAL  (A  -  H.  O  -  Z) 

PARAMETER  (NPT  =  150,  NPTT  =  50) 

COMMON  /NAME/  PRESIN.  TRAVIN.  ELECIN.  THERIN,  INTTIN 
COMMON  /DATA/  PRTIM  (1:NPT).  PRES  (1:NPT). 

1  TRTIM  (1:NPT).  TRAV  (1:NPT). 

2  VEL  (1:NPT).  ELTIM  (1:NPT),  ELEC  (1:NPT) 

COMMON  /OUT/  WFTIM  (NPT).  WFCUM  (NPT) 

COMMON  /NUMBER/  NPR.  NTR.  NEL.  NTH.  ITRVEL.  NWF 
COMMON  /START/  CHVOL.  AREA.  WFMASI.  RHOI. 

1  PROJI.  PROJW.  TOL.  TIMINC. 

2  TIMEND.  WAIR.  BM.  DERBM 

COMMON  /CURRENT/  TIME.  PRESCH.  TRAVPR. 

1  ELECEN.  TEMP.  RIMPET,  GAMMA, 

2  CO  VOL.  WFl.  WF2.  ELDEN,  SPMPR,  VELP,  BETA,  BASEPR,  CM, 

3  Z.  RHO.ENCOMP 

COMMON  /OPTIONS/  IBETA,  ICOM 

COMMON  /THERMO/  ELDENl  (LNPTT).  TEMPI  (LNPTT).  RIMPl  (LNPTT), 

1  G/VMMAl  (LNPTT).  COVOLl  (LNPTT) 

C 

CHARACTER*20  ETIN.  PRESIN.  TRAVIN.  ELEQN,  THERIN.  INITIN 
C 

WRITE  (*,  6000) 

6000  FORMAT(10X,’  Enter  the  file  containing  input  file  names.’) 

READ  (*.  5000)  ETIN 
OPEN  (15,  FILE  =  ETIN) 

REWIND  (15) 

C 

C  LINE  1:  PRESSURE 
C  LINE  2:  PROJECTILE  TRAVEL 
C  LINE  3:  ELECTRICAL  ENERGY  INPUT 

C  LINE  4:  THERMOCHEMICAL  VALUES  BASED  ON  ELECTRICAL  ENERGY  INPUT 
C  LINE  5:  INITIAL  DATA 

C  Each  file  contains  descriptor  on  line  1,  number  of  points  on  line  2, 

C  time  vs  parameter  on  other  lines  or  thermochems. 

READ  (15,  5000)  PRESIN 
READ  (15.  5000)  TRAVIN 
READ  (15.  5000)  ELEQN 
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READ  (15.  5(XX))  THERIN 
READ  (15.  5000)  INITIN 
CALL  OUTPUT  (1) 

CALL  PRESS  (1) 

CALL  TRAVEL  (1) 

CALL  ELECm  (1) 

CALL  INITIAL 

CALL  THERM  (1.  IREAC) 

CALL  OUTPUT  (2) 

C 

C  INITIALIZE  AT  TIME=0. 

C  SET  UP  ARRAY  FOR  WF  TO  FIND  DERIVATIVE 
TIME  =  0.0 
WF2  =  0. 

WFTLM  (1)  =  0.0 
WFCUM  (1)  =  0.0 
IEND  =  0 
NWF  =  1 
CALL  PRESS  (2) 

CALL  TRAVEL  (2) 

CALL  THERM  (2.  IREAC) 

CALL  OUTPUT  (3) 

C 

C  TAKE  FIRST  TIME  STEP 
TIME  =  TIME  +  TIMINC 
C 

C  INITIAL  GUESS  FOR  MASS  OF  FLUID  CONSUMED  IN  INITIAL 
C  ITERATE  TO  MATCH  THERMOCHEMS  WITH  MASS  CONSUMED 
1000  CONTINUE 
CALL  PRESS  (2) 

CALL  TRAVEL  (2) 

CALL  ELECTR  (2) 

M  =  0 
C 

C  IF  IREAC=0.  NO  REACTION  (EL  EN  DENSITY  NOT  LARGE  ENOUGH) 
C  IF  IREAC=1.  REACTION 
1010  CONTINUE 

CALL  THERM  (2.  IREAC) 

IF  (IREAC  .EQ.  0)  THEN 
CALL  SPMEAN 
CALL  OUTPUT  (3) 

TIME  =  TIME  +  TIMINC 
GO  TO  1000 
END  IF 
C 

CALL  STATE 
C 

C  M  IS  COUNTER-CONVERGE  BEFORE  100  STEPS  OR  SKIP 
IF  (M  .GT.  100)  THEN 
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WRITE  (*.  6010)  TIME 

6010  F0RMAT(1X.’  TIME='J^10.5,’  NO  CONVERGENCE’) 

TIME  =  TIME  +  TIMING 
GO  TO  1000 
END  IF 

IF  ((ABS  (WFl  -  WF2))  .LE.  TOL)  THEN 
C  FILL  ARRAY  FOR  WF  TO  FIND  DERIVATIVE 
NWF  =  NWF  +  1 
WFTIM  (NWF)  =  TIME 
WFCUM  (NWF)  =  WF2 
CALL  OUTPUT  (3) 

TIME  =  TIME  +  TIMING 

END  CONDITIONS 

IF  (TIME  .GT.  TIMEND)  GO  TO  1020 
IF  (WF2  .LE.  0.)  GO  TO  1020 
END  IF  ALL  FLUID  CONSUMED 

IF  ((WF2  .GE.  WFMASI)  .AND.  (lEND  .EQ.  0))  THEN 
WRITE  (*,  *)  ’  All  fluid  consumed.’ 

WRITE  (*,  *)  ’  Solution  continues  to  balance  energy.’ 

IEND=  1 
END  IF 

READY  TO  TAKE  ANO-PHER  STEP 
WFl  =  WF2 
GO  TO  1000 
ELSE 

WFl  =  WF2 
M  =  M+  1 
GO  TO  1010 
END  IF 

OOO  FORMAT  (A20) 

POST  PROCESSING-TAKE  DERIVATIVES 
1020  CONTINUE 
CALL  OUTPUT  (4) 

(XOSE  (26,  STATUS  =  ’KEEP’) 

CLOSE  (27.  STATUS  =  ’KEEP’) 

CLOSE  (28.  STATUS  =  ’KEEP’) 

CLOSE  (29.  STATUS  =  ’KEEP’) 

STOP 

END 

SUBROUTINE  PRESS  COPT) 

C 

IMPLICIT  REAL  (A  -  H,  O  -  Z) 

PARAMETER  (NPT  =  150.  NPTT  =  50) 

COMMON  /NAME/  PRESIN,  TRAVIN.  ELECIN.  THERIN,  INTTIN 
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COMMON  /DATA/  PRTIM  (1:NPT),  PRES  {1:NPT). 

1  TRTIM  (1  :NPT).  TRAV  ( 1  :NPT), 

2  VEL(lJSrPD.  ELTIM(1:NPT).  ELEC(1:NPT) 

COMMON  /NUMBER/  NPR,  NTR,  NEL,  NTH.  ITRVEL,  NWF 
COMMON  /CURRENT/  TIME,  PRESCH,  TRAVPR, 

1  ELECEN,  TEMP.  RIMPET.  GAMMA. 

2  COVOL.  WFl.  WF2.  ELDEN.  SPMPR.  VELP.  BETA.  BASEPR.  CM. 

3  Z.RHO.ENCOMP 
C 

CHARACTER*20  PRESIN.  TRAVIN.  ELEQN.  THERIN.  INITIN,  IDENT 
C 
C 

C  IOPT=l  READ  RLE 
C  TIME  IN  MS.  PRESSURE  IN  MPA 
C  CHANGE  TIME  TO  SEC 

C  CHANGE  PRESSURE  FROM  MPA  TO  CONSISTENT  UNTFS  WITH  CGS  SYSTEM 
C  BY  MULTIPLYING  BY  L0E7  CONVERSION  CONSTANT 
C 

IF  aOPT  .EQ.  2)  GO  TO  1000 

OPEN  (21.  FILE  =  PRESIN.  STATUS  =  ‘OLD*.  lOSTAT  =  lOS,  ERR  =  1010) 
REWIND  (21) 

1010  CONTINUE 

IF  aOS  .NE.  0)  WRITE  (*.  6000)  lOS 
6000  FORMAT(’  ERROR  OPENING  PRESSURE  FILE: ’.15) 

READ  (21.  5000)  IDENT 
WRITE  (*.  6010)  PRESIN 
C  WRITE(*.620)  IDENT 
READ  (21.  *)  NPR 
DO  1020  I  =  1.  NPR 
READ  (21.  *)  PRTIM  G),  PRES  (I) 

C  CHANGE  TO  SEC  AND  CGS  SYSTEM 
PRTIM  (I)  =  PRTIM  (I)*1.0E-3 
PRES  (I)  =  PRES  a)*L0E7 
1020  CONTINUE 
GO  TO  1030 
C 

C  IOPT=2  INTERPOLATE 
1000  CONTINUE 

TIMTOL  =  ABS  (TIME  -  .lE-10) 

IF  ((TIME  .EQ.  0.0)  .or.  (TIMTOL  .EQ.  0.))  THEN 
PRESCH  =  PRES  (1) 

CALL  SPMEAN 
END  IF 

DO  1040  I  =  2.  NPR 

IF  (TIME  .EQ.  PRTIM  H))  PRESCH  =  PRES  (I) 

IF  (TIME  .GT.  PRTIM  G  -  1))  THEN 

pREsai  =  PRES  a  - 1)  +  (PRES  (I)  -  PRES,  a  -  1))* 

1  ((TIME  -  PRTIM  (I  -  1))/(PRTIM  (I)  -  PRTIM  Q.  -  1))) 

END  IF 
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1040  CONTINUE 
C 

5000  FORMAT  (A80) 

6010  FORMATC  Pressure-time  file;*.A20) 

8000  FORMATC  Pressure-time  ident:’,A20/) 

1030  CONTINUE 
RETURN 
END 

0  «4i*4i************4i4i*****4i*««4i4i**** 

SUBROUTINE  TRAVEL  OOPT) 

C 

C  DETERMINE  PROJECTILE  TRAVEL 

C  CALL  SUBROUTINE  TO  DETERMINE  PROJECTILE  VELOCITY 
IMPLICIT  REAL  (A  -  H.  O  -  Z) 

PARAMETER  (NPT  =  150.  NPTT  =  50) 

COMMON  /NAME/  PRESIN.  TRAVIN.  ELECIN.  THERIN.  INTTIN 
COMMON  /DATA/  PRTIM  (1:NPT).  PRES  (1:NPT). 

1  TRTIM  (1:NPT),  TRAV  (1:NPT). 

2  VEL(1:NPT).  ELTIM(1:NPT).  ELEC(1:NPT) 

COMMON  /NUMBER/  NPR.  NTR.  NEL.  NTH.  ITRVEL.  NWF 
COMMON  /CURRENT/  TIME,  PRESCH.  TRAVPR. 

1  ELECEN.  TEMP.  RIMPET.  GAMMA. 

2  CO  VOL.  WFl,  WF2.  ELDEN.  SPMPR.  VELP,  BETA.  BASEPR,  CM. 

3  Z.  RHO.  ENCOMP 

COMMON  /OPTIONS/  IBETA,  ICOM 
C 

CHARACTER*20  PRESIN.  TRAVIN.  ELEQN.  THERIN,  INITIN,  IDENT 
C 

IF  (GOPT  .NE.  1)  .AND.  (lOPT  .NE.  2)) 

1  WRITE  (*,  *)  ’  Error  in  Subroutine  Travel.’ 

C 

C  IOPT=l  READ  HLE 

C  IF  2ND  LINE  IS  1:  TIME  IN  MS.  TRAVEL  IN  CM 

C  IF  2ND  LINE  IS  2:  TIME  IN  MS.  TRAVEL  IS  IN  CM,  VELOCITY  IN  CM/S 
C  CHANGE  TIME  TO  SEC 
C  TRAVEL  MUST  BE  IN  CM 
C  VELOCITY  IN  CM/S 
C 

IF  GOPT  .EQ.  2)  GO  TO  1000 

OPEN  (21,  FILE  =  TRAVIN,  STATUS  =  ’OLD’,  lOSTAT  =  lOS,  ERR  =  1010) 
REWIND  (21) 

1010  CONTINUE 

IF  GOS  .NE.  0)  WRITE  (•.  6000)  lOS 
6000  FORMAT(’  ERROR  OPENING  TRAVEL  FILE:’.I5) 

READ  (21.  5000)  IDENT 
WRITE  (*.  6010)  TRAVIN 
C  WRITE(*.620)  IDENT 
C  ITRVEL=1  FOR  TIME,  TRAVEL  PAIRS 
C  ITRVEL=2  FOR  TIME,  TRAVEL.  VELOQTY  PAIRS 
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READ  (21.  *)  ITRVEL 

IF  (ITRVEL  .EQ.  1)  WRITE  (*.  6020) 

IF  (ITRVEL  .EQ.  2)  WRITE  (*.  6030) 

IF  ((ITRVEL  .NE.  1)  .AND.  (ITRVEL  .NE.  2))  THEN 
WRITE  (*.  6040) 

STOP 
END  IF 

READ  (21.  *)  NTR 
DO  1020  I  =  1.  NTR 

READ  TIME-TRAVEL  IF  ITRVEL=1 
CHANGE  TIME  TO  SEC 
**IBHVG2 

CHANGE  TRAVEL  TO  CM  (GIVEN  IN  .1  MM).  VEL  TO  CM/SEC 
IF  (ITRVEL  .EQ.  1)  THEN 
READ(21.*)  TRTIMa).TRAVa) 

TRTIMa)=TRTIM(I)*  1  .OE-3 
TRAV(I)=TRAVa)*0.01 
ENDIF 

READ  TIME-TRAVEL- VELOCITY  IF  ITRVEL=2 
IF  (ITRVEL  .EQ.  2)  THEN 
READ(21.*)  TRTIMa).TRAV(I).VEL(I) 
TRTIMa)=TRTIM(I)*  1  .OE-3 
TRAV(I)=TRAVa)*0.01 
VELa)=VEL(I)*0.1 
ENDIF 

TIME  IN  MS.  TRAVEL  GIVEN  IN  CM 
IF  (ITRVEL  .EQ.  1)  THEN 
READ  (21.  •)  TRTIM  (I).  TRAV  © 

TRTIM  0)  =  TRTIM  a)*1.0E-3 
END  IF 

020  CONTINUE 
GO  TO  1030 

IOPT=2  INTERPOLATE 
000  CONTINUE 

IF  (TIME  .EQ.  0.)  THEN 
TRAVPR  =  TRAV  (1) 

VELP  =  0.0 
ENDIF 


ITRVEL=1  HAVE  JUST  TIME-TRAVEL 
IF  (ITRVEL  .EQ.  1)  THEN 
DO  1040  I  =  2.  NTR 
IF  (TIME  .EQ.  TRTIM  G))  THEN 
TRAVPR  =  TRAV  0) 
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CALL  VELPR  (I) 

GO  TO  1030 
END  IF 

IF  (TIME  .GT.  TRTIM  G  -  1))  THEN 
TRAVPR  =  TRAV  0  -  1)  +  (TRAV  (I)  -  TRAV  0  -  1))* 

1  ((TIME  -  TRTIM  G  -  1))/(TRTIM  G)  -  TRTIM  G  *  1))) 

CALL  VELPR  (I) 

END  IF 

1040  CONTINUE 
END  IF 
C 

C  ITRVEL=2  HAVE  TIME-TRAVEL- VELOCITY 
IF  GTRVEL  .EQ.  2)  THEN 
DO  1050  I  =  2.  NTR 
IF  (TIME  .EQ.  TRTIM  G))  THEN 
TRAVPR  =  TRAV  G) 

VELP  =  VEL  (I) 

GO  TO  1030 
END  IF 

IF  (TIME  .GT.  TRTIM  G  *  1))  THEN 
TRAVPR  =  TRAV  G  -  1)  +  (TRAV  (I)  -  TRAV  G  -  1))* 

1  ((TIME  -  TRTIM  G  -  1))/(TRTIM  G)  -  TRTIM  G  -  1))) 

VELP  =  VEL  (I  •  1)  +  (VEL  G)  -  VEL  G  *  1))* 

1  ((TIME  -  TRTIM  G  -  1))/(TRTIM  G)  -  TRTIM  G  -  1))) 

END  IF 

1050  CONTINUE 
END  IF 
C 

5000  FORMAT  (A80) 

6010  FORMATC  Proj  motion-time  file:’^0) 

80(X)  FORMATC  Proj  motion-time  ident:’,A20y) 

6020  FORMATC  Proj  travel-time  flle.’y) 

6030  FORMAT(’  Proj  velocity-time  file.’/) 

6040  FORMAT(’  Error-unknown  type  file.’/) 

1030  CONTINUE 
RETURN 
END 

SUBROUTINE  ELECTR  GOPT) 

C 

IMPLICIT  REAL  (A  -  H,  O  -  Z) 

PARAMETER  G^FT  =  150,  NPTT  =  50) 

COMMON  /NAME/  PRESIN,  TRAVIN,  ELECIN.  THERIN,  INITIN 
COMMON  /DATA/  PRTIM  (1:NPT).  PRES  (1;NPT). 

1  TRTIM  (1  :NPT),  TRAV  (1 

2  VEL  (IJ^.  ELTIM  (1:NPT).  ELEC  (1:NPT) 

COMMON  /NUMBER/  NPR,  NTR.  NEL,  NTH.  ITRVEL.  NWF 
COMMON  /CURRENT/  TIME.  PRESCH.  TRAVPR, 

1  ELECEN,  TEMP.  RIMPET,  G/VMMA, 
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2  COVOL.  WFl,  WF2.  ELDEN,  SPMPR,  VELP,  BETA.  BASEPR,  CM. 

3  Z.  RHO,  ENCOMP 

COMMON  /OPTIONS/  IBETA,  ICOM 
C 

CHARACTER*20  PRESIN,  TRAVIN,  ELEQN.  THERIN,  INITIN,  IDENT 
C 

IF  (GOPT  .NE.  I)  .AND.  (lOPT  .NE.  2)) 

1  WRITE  (*,  *)  ’  Error  in  Subroutine  Electr.’ 

C 

C  IOPT=I  READ  HLE 
C  ELECTRICAL  ENERGY  IS  CUMULATIVE 
C  TIME  IN  MS.  ELECTRICAL  ENERGY  IN  kJ 
C  CHANGE  TIME  TO  SEC 
C 

IF  (lOPT  .EQ.  2)  GO  TO  1000 

OPEN  (21,  FILE  =  ELECIN.  STATUS  =  ’OLD*.  lOSTAT  =  lOS.  ERR  =  1010) 
REWIND  (21) 

1010  CONTINUE 

IF  (lOS  .NE.  0)  WRITE  (*.  6000)  lOS 
6000  FORMATC  ERROR  OPENING  ELECTRICAL  ENERGY  RLE: *.15) 

READ  (21.  5000)  IDENT 
WRITE  (♦.  6010)  ELECIN 
C  WRITE(*.620)  IDENT 
READ  (21.  *)  NEL 
DO  1020  I  =  1.  NEL 
READ  (21.  *)  ELTIM  G),  ELEC  (I) 

ELTIM  G)  =  ELTIM  (I)*1.0E-3 
1020  CONTINUE 
GO  TO  1030 
C 

C  IOPT=2  INTERPOLATE 
1000  CONTINUE 

IF  (TIME  .EQ.  0.)  ELECEN  =  ELEC  (1) 

DO  1040  I  =  2.  NEL 

IF  (TIME  EQ.  ELTIM  G))  ELECEN  =  ELEC  G) 

IF  (TIME  .GT.  ELTIM  G  -  1))  THEN 
ELECEN  =  ELEC  G  -  1)  +  (ELEC  G)  -  ELEC  (I  -  1))* 

1  ((TIME  -  ELTIM  (I  -  1))/(ELTIM  (I)  -  ELTIM  G  -  1))) 

END  IF 

1040  CONTINUE 
C 

5000  FORMAT  (A80) 

6010  FORMATC  Electrical  energy-time  file:’.A20) 

SCXX)  FORMAT(’  Electrical  energy-time  ident:’,A20y) 

1030  CONTINUE 
RETURN 
END 

SUBROUTINE  INITIAL 
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IMPUCrr  REAL  (A  .  H.  O  -  Z) 

PARAMETER  (MPT  =  150.  NPTT  =  50) 

COMMON  /NAME/  PRESIN.  TRAVIN.  ELECIN.  THERIN.  INITIN 
COMMON  /START/  CHVOL.  AREA.  WFMASI.  RHOI. 

1  PROJI.  PROJW,  TOL.  TIMINC. 

2  TIMEND.  WAIR,  BM.  DERBM 

COMMON  /CURRENT/  TIME.  PRESCH.  TRAVPR. 

1  ELECEN.  TEMP.  RIMPET,  GAMMA. 

2  COVOL.  WFl.  WF2.  ELDEN.  SPMPR.  VELP.  BETA.  BASEPR.  CM. 

3  Z,RHO.ENCOMP 

COMMON  /OPTIONS/  IBETA.  ICOM 
C 

CHARACTER*20  PRESIN.  TRAVIN.  ELEQN.  THERIN.  INITIN.  IDENT 
C 

OPEN  (21.  FILE  =  INITIN.  STATUS  =  'OLD'.  lOSTAT  =  lOS.  ERR  =  1000) 
REWIND  (21) 

1000  CONTINUE 

IF  (lOS  .NE.  0)  WRITE  (*.  6000)  lOS 
6000  FORMATC  ERROR  OPENING  INITIAL  DATA  FILE: ’.15) 

C 

C  CHAMBER  VOLUME  IN  CM'^S 

C  INITIAL  PROJECTILE  POSITION  IN  CM  WITH  0.0  AT  ORIGIN  OF  TUBE 
C  TUBE  DIAMETER  IN  CM-AREA  IN  CM'^2 
C  PROJECTILE  MASS  IN  KG 
C  INITIAL  MASS  OF  WORKING  FLUID  IN  GRAMS 
DENSITY  OF  WORKING  FLUID  IN  G/CM''3 
TIME  INCREMENT  IN  SEC 
READ  (21.  5000)  IDENT 
WRITE  (•,  6010)  INITIN 
C  WRITE(*.620)  IDENT 
READ  (21.  *)  CHVOL 
WRITE  (*.  6020)  CHVOL 
C  TUBE  DIAMETER 
READ  (21.  •)  DIAM 

AREA  =  3.141592654*(DIAM/2.)*(DIAM/2.) 

WRITE  (•.  6030)  DIAM.  AREA 

C  PROJECTILE  INITIAL  POSITION  AND  MASS  (CHANGE  KG  TO  G) 

READ  (21.  •)  PROJI 

IF  (PROJI  .NE.  0.0)  WRITE  (*.  •) 

1  '  Change  needed-proj  not  zeroed.' 

READ  (21.  •)  PROJW 
WRITE  (•.  6040)  PROJI.  PROJW 
PROJW  «  PROJW*1000. 

C  WORKING  FLUID  INITIAL  MASS 
READ  (21.  •)  WFMASI 
WRITE  (*.  6050)  WFMASI 
READ  (21.  •)  RHOI 
WRITE  (•.  6060)  RHOI 
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RHO  =  RHOI 

C  TOLERANCE  IN  ITERATION 
READ  (21.  •)  TOL 
WRITE  (*.  6070)  TOL 
C  TIME  IN  SEC 

READ  (21.  •)  TIME 
READ  (21.  •)TIMINC 
READ  (21.  *)TIMEND 
WRITE  (*.  6080)  TIME.  TIMINC.  TIMEND 
C  INITIAL  GUESS  FOR  MASS  OF  WORKING  FLUID  CONSUMED 
READ  (21.  •)  WFl 
WRITE  (*.  6090)  WFl 
C  READ  INITIAL  VALUE  OF  BETA 
C  OPTIONS  FOR  BETA: 

C  BETA  IS  INITIAL  VALUE  OF  WORKING  FLUID  MOVING-- 
C  LAGRANGE  ASSUMPTION  IS  BETA=1.  I.E.  ALL  FLUID  MOVES 

C  IF  BETA=0..  FLUID  IS  AT  REST 

C  CHANGE  FRACTION  OF  BETA  MOVING  IF  IBETA=1 
C  IF  IBETA=0.  LET  BETA=CONSTANT  (INITIAL  VALUE) 

C  IF  IBETA=1.  LET  BETA=Z  (FRACTTION  OF  WF  REACTED) 

READ  (21.  *)  BETA.  IBETA 
WRITE  (*.  6100)  BETA 
IF  (IBETA  .EQ.  1)  THEN 
WRITE  (•.  6110)  IBETA 
Z  =  0. 

END  IF 

C  MASS  OF  AIR  INITIALLY  IN  CHAMBER 
READ  (21.  *)  WAIR 
WRITE  (*.  6120)  WAIR 
C 

C  OPTION  -  COMPRESSION  OF  WORKING  FLUID 
C  IF  ICOM=0.  DO  NOT  CONSIDER  COMPRESSION 
C  IF  IC0M=1.  CONSIDER  COMPRESSION 
READ  (21.  •)  ICOM 

IF  (OCOM  .NE.  0)  .AND.  (ICOM  .NE.  1))  WRITE  (♦.  *)  ’Error-ICOM’ 

C  BM=BULK  MODULUS  OF  FLUID  AT  ONE  ATMOSPHERE  (MPa) 

C  DERBM=DERIVATIVE  OF  BULK  MODULUS  WITH  RESPECT  TO  PRESSURE 
IF  OCOM  £Q.  1)  THEN 
READ  (21.  *)  BM.  DERBM 
WRITE  (*.  6130)  BM.  DERBM 
C  CHANGE  TO  CGS 
BM  =  BM*1.0E7 
END  IF 
C 

5000  FORMAT  (A80) 

6010  FORMATC  IniUal  data  file:’.A20) 

8000  FORMATC  IniUal  data  ident:’,A20) 

6020  FORMAT(’  (Chamber  volume  (cm''3)=’.F15.5) 

6030  FORMAT(’  Tube  diameter  (cm)=’JF6.I.2x.’  Area  (cm^2)=’J^12.5) 
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6040  FORMATC  Projectile  position  (cm)=*f  10.5, 

1  2XJ'  Projectile  mass  (kg)=’ J^IO.2) 

6050  FORMATC  Working  Huid  mass  (g)=T10J) 

6060  FORMATC  Initial  working  fluid  density  (g/cm^3)=’,F10.5) 

6070  FORMATC  Tolerance  required  in  mass  consumed=*,F10.5) 

6080  FORMATC  'Hnie  start  (sec)=’ J^lO.Sy’  Time  increment  (sec)=’, 

1  F10.8y’  Time  end  (sec)=’J^10.8) 

6090  FORMAT(*  Initial  guess  for  mass  consumed  (g)=’,F10.5) 

6100  FORMATC  Initial  value  of  beta,  fraction  of  wf  moving  (-)=’ J^IO.5) 

6110  FORMATC  Option:  ibeta=‘J5y. 

1 '  If  ibeta=0  fraction  moving  is  constantV. 

2’  If  ibeta=l.  fraction  moving=fraction  decomposed') 

6120  FORMATC  Mass  of  air  initiaUy  in  chamber  (g)=’J^10.5) 

6130  FORMATC  Bulk  modulus  of  WF  at  one  atmosphere  (MPa)=’JF10.5y, 
r  Derivative  of  bulk  modulus  wrt  pressure  (->=’,F10.5y) 

C 

CXOSE  (21) 

RETURN 

END 

Q  4i4i*4i*«i4i***4i*4i#4(4i*4i*******«r«4i*4i*«* 

SUBROUTINE  THERM  GOPT.  IREAQ 
C 

IMPLICIT  REAL  (A  -  H.  O  -  Z) 

PARAMETER  (NPT  =  150,  NPTT  =  50) 

COMMON  /NAME/  PRESIN,  TTIAVIN.  ELECIN,  THERIN,  INITIN 
COMMON  /NUMBER/  NPR,  NTR,  NEL,  NTH,  ITRVEL.  NWF 
COMMON  /START/  CHVOL,  AREA,  WFMASI.  RHOI, 

1  PROJI,  PROJW,  TOL,  TIMINC, 

2  TIMEND.  WAIR,  BM.  DERBM 

COMMON  /CURRENT/  TIME.  PRESCH,  TRAVPR, 

1  ELECEN,  TEMP.  RIMPET,  GAMMA, 

2  COVOL.  WFl,  WF2,  ELDEN,  SPMPR,  VELP,  BETA.  BASEPR,  CM. 

3  Z.  RHO,  ENCOMP 

COMMON  /OPTIONS/  IBETA,  ICOM 

COMMON /THERMO/ ELDENl  (1:NPTT).  TEMPI  (1:NPTT),  RIMPl  (1:NPTT), 
1  GAMMAl  (LNPTT),  COVOLl  (UNFIT) 

C 

CHARACTER*20  PRESIN,  TRAVIN,  ELEQN.  THERIN.  INITIN,  IDENT 
C 

IF  (GOPT  .NE.  1)  .AND.  (lOPT  .NE.  2)) 

1  WRITE  (*,  *)  ’  Error  in  Subroutine  Electr.’ 

C 

C  IOPT=l  READ  HLE 

IF  GOPT  .EQ.  2)  GO  TO  1000 

C  BASED  ON  AMOUNT  OF  ELECTRICAL  ENERGY  INPUT 

OPEN  (21,  FILE  =  THERIN,  STATUS  =  'OLD'.  lOSTAT  =  lOS,  ERR  =  1010) 
REWIND  (21) 

1010  CONTINUE 

IF  GOS  .NE.  0)  WRITE  (•,  6000)  lOS 
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6000  FORMATC  ERROR  OPENING  THERMOCHEM  DATA  FILE:  *.15) 

C 

C  ELECTRICAL  ENERGY  DENSITY  (Elec  enAnass  fluid  in  chamber)  IN  kJ/G 
C  TEMPERATURE  IN  K 
C  IMPETUS  IN  J/G 

C  GAMMA,  RATIO  OF  SPECIHIC  HEAT.  UNITLESS 
C  COVOLUME  IN  CM'^3/G 

C  CONVERT  JOULES  TO  CGS:  1  JOULE=1.0E7  G-CM''2/S'^2 
READ  (21.  5000)  IDENT 
WRITE  (*.  6010)  THERIN 
C  WRITE(*.620)  IDENT 
READ  (21.  *)  NTH 
DO  1020  I  =  1.  NTH 

READ  (21,  *)  ELDENl  G).  TEMPI  G).  RIMPl  G). 

1  GAMMAl  G).  COVOLl  (1) 

C  ELECTRICAL  ENERGY  DENSITY  IN  kJ/G 
RIMPl  G)  =  RIMPl  G)*1.0E+7 
1020  CONTINUE 
GO  TO  1030 
C 

1000  CONTINUE 

C  ELECTRICAL  ENERGY  DENSITY=CUMULATIVE  ELEC  ENERGY/CUMULATIVE  WF  REACTED 
IF  (WFl  .LT.  0.0)  THEN 
WRITE  (*.  6020) 

END  IF 

C  ASSUME  FIRST  LINE  OF  TABLE  IS  0. 

IF  ((WFl  .LE.  0.)  .OR.  (TIME  .EQ.  0.))  THEN 
TEMP  =  TEMPI  (1) 

RIMPET  =  RIMPl  (1) 

GAMMA  =  GAMMAl  (1) 

COVOL  =  COVOLl  (1) 

GO  TO  1030 
END  IF 
C 

IF  (WFl  .GT.  0.)  ELDEN  =  ELECEN/WFl 
C 

C  IF  WF  CONSUMED.  FIX  THERMOCHEMS  AT  FINAL  VALUE 
IF  (WFl  .GT.  WFMASI)  ELDEN  =  ELECEN/WFMASI 
C 

C  IF  EL  EN  DENSITY  NOT  LARGE  ENOUGH  FOR  REACTION, 

C  COMMUNICATE  TO  MAIN  WITH  IREAC=0  AND  SKIP  THIS  STEP 
IF  (ELDEN  .LT.  ELDENl  (1))  THEN 
IREAC  =  0 
GO  TO  1030 
END  IF 
C 

DO  1040  I  =  2.  NTH 
IF  (ELDEN  .EQ.  ELDENl  G))  THEN 
TEMP  =  TEMPI  (I) 
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RIMPET  =  RIMPl  a) 

GAMMA  =  GAMMAl  (I) 

COVOL  =  COVOLl  (I) 

END  IF 

IF  (ELDEN  .GT.  ELDENl  G  -  D)  THEN 
TEMP  =  TEMPI  a  '  1)  +  (TEMPI  (1)  -  TEMPI  0  -  D)* 

1  ((ELDEN  -  ELDENl  Q  -  1))/(ELDEN1  (I)  -  ELDENl  G  -  1))) 
RIMPET  =  RIMPl  G  -  1)  +  GGMPl  G)  -  RIMPl  G  -  D)* 

1  ((ELDEN  -  ELDENl  G  *  1))/(ELDEN1  G)  -  ELDENl  G  -  1))) 

GAMMA  =  GAMMAl  G  -  1)  +  (GAMMAl  (I)  -  GAMMAl  G  -  1))* 
1  ((ELDEN  -  ELDENl  G  -  1))/(ELDEN1  G)  -  ELDENl  (I  -  1))) 

COVOL  =  COVOLl  (I  -  1)  +  (COVOLl  G)  -  COVOLl  (1  -  1))* 

1  ((ELDEN  .  ELDENl  G  -  1))/(ELDEN1  G)  -  ELDENl  G  -  D)) 
END  IF 
IREAC=  1 
1040  CONTINUE 
C 

C  IF  FALL  OUTSIDE  ARRAY 

IF  (ELDEN  .GT.  ELDENl  (NTH))  THEN 
TEMP  =  TEMPI  (NTH) 

RIMPET  =  RIMPl  (NTH) 

GAMMA  =  GAMMAl  (NTH) 

COVOL  =  COVOLl  (NTH) 

WRITE  (26,  *)  ’  OUTSIDE  THERM  ARRAY’ 

END  IF 
C 

5000  FORMAT  (A80) 

6010  FORMAT(’  Thermochemical  data  file:’,A20) 

80(X)  FORMATC  Thermochemical  data  ident:’,A20y) 

6020  FORMATC  Error  in  Subroutine  Electr— El  En  Density  negative.’) 

C 

1030  CONTINUE 
RETURN 
END 

^  **4i4i4r**iii***4i**4r**4>****4i**4,4,*4,4,4,4,« 

SUBROUTINE  STATE 
C 

IMPLICIT  REAL  (A  -  H.  0  -  Z) 

COMMON  /START/  CHVOL,  AREA,  WFMASI,  RHOI, 

1  PROJI,  PROJW,  TOL,  TIMINC, 

2  TIMEND,  WAIR,  BM,  DERBM 

COMMON  /CURRENT/  TIME,  PRESCH,  TRAVPR, 

1  ELECEN,  TEMP,  RIMPET,  GAMMA, 

2  COVOL,  WFl,  WF2,  ELDEN,  SPMPR,  VELP,  BETA,  BASEPR,  CM, 

3  Z,RHO,ENCOMP 
COMMON  /OPTIONS/  IBETA,  ICOM 

C 

C  DETERMINE  SPACE-MEAN  PRESSURE 
CALL  SPMEAN 


n  n 


c 

C  IF  ICOM=0.  DO  NOT  CONSIDER  COMPRESSION  OF  WF 
C  IF  ICOM=l,  CONSIDER  COMPRESSION  OF  WF 

C  ENCOMP  IS  ENERGY  PER  UNIT  MASS  STORED  BY  COMPRESSED  FLUID 
C  -NEED  TO  MULTIPLY  BY  AMT  OF  FLUID  UNV.M>ORIZED 
IF  aCOM  £Q.  0)  THEN 
RHO  =  RHOI 
ENCOMP  =  0.0 
END  IF 

IF  aCOM  £Q.  1)  CALL  COMPRES 

CUMULATIVE  FLUID  CONSUMED  FROM  ENERGY  EQN,  P  105,  BK  9 
GAMl  =  GAMMA  -  1. 

DIV  =  RIMPET/GAMl  +  SPMPR*(COVOL/GAMl  -  l./(RHO*GAMl)) 

C 

TOP  =  (CHVOL  -  WFMASI/RHO  +  TRAVPR*AREA)*(SPMPRyGAMl) 

1  +  VELP*VELP*((U6.)*CM*PROJW  +  .5*PROJW) 

2  +  (1.  -  Z)*WFMASI*ENCOMP 
C 

Z  =  (L/WFMASI)*(TOP/DIV) 

WF2  =  Z*WFMASI 
C 

RETURN 

END 

SUBROUTINE  SPMEAN 
C 

IMPLICIT  REAL  (A  -  H.  O  -  Z) 

COMMON  /START/  CHVOL,  AREA.  WFMASI,  RHOI. 

1  PROJI,  PROJW,  TOL,  TIMINC. 

2  TIMEND,  WAIR,  BM.  DERBM 

COMMON  /CURRENT/  TIME,  PRESCH.  TRAVPR, 

1  ELECEN,  TEMP.  RIMPET,  GAMMA, 

2  COVOL,  WFl,  WF2,  ELDEN,  SPMPR,  VELP,  BETA,  BASEPR,  CM. 

3  Z,  RHO,  ENCOMP 

COMMON  /OPTIONS/  IBETA,  ICOM 
C 

C  SPACE  MEAN  PRESSURE  EQUALS  CHAMBER  PRESSURE  (IF  EXPERIMENTAL) 
C  SPMPR=PRESCH 
C 

C  SPACE-MEAN  PRESSURE  FROM  LAGRANGE  GRADIENT 
C  ASSUME  PRESCH  IS  A  PRESSURE  MEASUREMENT  AT  THE  BREECH 
C  INCREASE  PROJECTILE  MASS  BY  5%  TO  ACCOUNT  FOR  FRICTION 
C  CHARGE  TO  MASS  RATIO,  CM.  IS  INITIAL  MASS  OF  WORKING  FLUID 
C  DIVIDED  BY  MASS  OF  PROJECTILE 
C  projw=1.05*PROJW 
C  For  test  case,  no  friction 
ALPHA  =  WF2/WFMASI 
C 
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C  CONSIDER  BETA=FRACnON  OF  FLUID  WHICH  MOVES 
C  IF  BETA=0.,  ALL  THE  FLUID  IS  AT  REST 
C  IF  BETA=L.  ALL  THE  WF  MOVES  (LAGRANGE  ASSUME) 

IF  (IBETA  .EQ.  1)  BETA  =  Z 

CM  =  (WFMASI*(ALPHA  +  BErA*(L  -  ALPHA))  +  WAIR)/PROJW 
IF  (CM  1£.  0.0)  WRITE  (*,  •)  ’  Error  in  Spmcan-C/M  <=  0.’ 

SPMPR  =  ((1.  +  iy3.*CM)/(L  +  172.*CM))*PRESCH 
BASEPR  =  (17(1.  +  173.*CM))*SPMPR 
C 

RETURN 

END 

SUBROUTINE  COMPRES 

DETERMINE  DENSITY  OF  COMPRESSED  WORKING  FLUID 
DETERMINE  ENERGY  OF  COMPRESSION 

IMPLICIT  REAL  (A  -  H.  O  -  Z) 

PARAMETER  (NPT  =  150.  NPTT  =  50) 

COMMON  /START/  CHVOL,  AREA,  WFMASI.  RHOI, 

1  PROJI,  PROJW.  TOL.  TIMINC, 

2  TIMEND,  WAIR.  BM.  DERBM 
COMMON  /CURRENT/  TIME,  PRESCH,  TRAVPR, 

1  ELECEN.  TEMP,  RIMPET,  GAMMA, 

2  COVOL,  WFl.  WF2.  ELDEN,  SPMPR.  VELP,  BETA,  BASEPR,  CM, 

3  Z.RHO.ENCOMP 
COMMON  /OPTIONS/  IBETA.  ICOM 
COMMON  /DATA/  PRTIM  (1:NPT),  PRES  (LNPT), 

1  TRTIM(l:NPT),TRAV(ld^. 

2  VEL(1:NPT).  ELTIM(1:NPT),  ELEC(1:NPT) 

LET  BULK  MODULUS  BE  AT  MOST  A  LINEAR  FCN  OF  PRESSURE 
IF  (BM  iE.  0.0)  WRITE  (*,  *)  ’ERROR-BM<=0.0’ 

IF  (DERBM  .EQ.  0.0)  THEN 
RHO  =  RHOI*(l.  +  (PRESCH  -  PRES  (1))/BM) 

ELSE 

RHO  =  RHOI*(L  +  (DERBM/BM)*(PRESCH  -  PRES  (!)))**(! 7DERBM) 
END  IF 

ENERGY  PER  UNIT  MASS  STORED  BY  COMPRESSED  FLUID 
IF  ((DERBM  £Q.  0.0)  .OR.  (DERBM  £Q.  1.0))  THEN 
ENCOMP  =  .  l.*(PRESCH  -  PRES  (l))/RHO  +  (BM/RHOI)* 

1  LOG  (1.  +  (PRESCH  -  PRES  (1))/BM) 

ELSE 

ENCOMP  =  -  l.*(PRESCH  -  PRES  (l))/RHO  +  (BM  + 

1  DERBM*(PRESCH  -  PRES  (1))) 

2  /((DERBM  -  l.)*RHO)  -  BM/((DERBM  -  l.)*RHOI) 

END  IF 


57 


uuu  uu  u  u  uuuu 


RETURN 

END 

Q  ^,*!^**m*****************i********* 

SUBROUTINE  VELPR  (I) 

DETERMINE  PROJECTILE  VELOCITY 

IMPLICIT  REAL  (A  -  H.  O  -  Z) 

PARAMETER  (NFT  =  150.  NPTT  =  50) 

COMMON  /DATA/  PRTIM  (1:NPT).  PRES  (1:NPT). 

1  TRTIM  (LNPT),  TRAV  (1 

2  VEL  (LNPT).  ELTIM  (LNPT).  ELEC  (LNPT) 

COMMON  /NUMBER/  NPR.  NTR.  NEL.  NTH.  ITRVEL.  NWF 
COMMON  /START/  CHVOL.  AREA.  WFMASI.  RHOI. 

1  PROJI.  PROJW.  TOL.  TIMINC. 

2  TIMEND.  WAIR.  BM.  DERBM 
COMMON  /CURRENT/  TIME.  PRESCH,  TRAVPR. 

1  ELECEN.  TEMP.  RIMPET.  GAMMA. 

2  CO  VOL.  WFl.  WF2.  ELDEN.  SPMPR.  VELP.  BETA.  BASEPR,  CM. 

3  Z.  RHO.  ENCOMP 
COMMON  /OPTIONS/  IBETA.  ICOM 

CENTERED  DIFFERENCE 
IF  a  EQ.  NTR)  THEN 
VELP  =  (TRAV  a)  -  TRAV  (I  -  2))/(2.* 

1  (TRTIM  a)  -  TRTIM  (I  -  1))) 

ELSE 

VELP  =  (TRAV  a  +  1)  -  TRAV  G  *  l))/(2.* 

1  (TRTIM  0)  -  TRTIM  (I  -  1))) 

END  IF 

RETURN 
END 

SUBROUTINE  DER  (TIME.  lOPT.  AMT.  AMTDER) 

FINDS  DERIVATIVES:  USED  FOR  dE/dt  and  dWF/dt 
ASSUMES  CONSTANT  STEP  SIZE  FOR  TIME 

IMPLICIT  REAL  (A  -  H.  O  -  Z) 

PARAMETER  (NPT  =  150.  NPTT  =  50) 

COMMON  /DATA/  PRTIM  (LNPT).  PRES  (LNPT). 

1  TRTIM  (LNPT).  TRAV  (10^. 

2  VEL  (LNPT).  ELTIM  (LNPT).  ELEC  (IJSrPT) 

COMMON  /OUT/  WFTIM  (NPT).  WFCUM  (NPT) 

COMMON  /NUMBER/  NPR.  NTR.  NEL.  NTH.  ITRVEL.  NWF 

C 

C  IF  IOPT=l.  FIND  DER  OF  ENERGY 
C  IF  IOPT=2.  FIND  DER  OF  WF 
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IF  aOPT  .EQ.  1)  GO  TO  1000 
IF  aOPT  .EQ.  2)  GO  TO  1010 
C 

IF  ((ABS  (TIME  -  l.E-10))  £Q.  0.0)  THEN 
AMT  =  0.0 
AMTDER  =  0.0 
GO  TO  1020 
END  IF 
C 

C  IOPT=l,  DERIVATIVE  OF  ELEC  EN 
1000  CONTINUE 

DO  1030  I  =  2,  NEL 
IF  (TIME  .EQ.  ELTIM  (I))  THEN 
AMT  =  ELEC  (I) 

IMARK  =  I 
GO  TO  1040 
END  IF 

IF  (TIME  .GT.  ELTIM  -  1))  THEN 
AMT  =  ELEC  (I  -  1)  +  (ELEC  H)  -  ELEC  G  -  1))* 

1  ((TIME  -  ELTIM  (I  -  1))/(ELTIM  (I)  -  ELTIM  G  *  1))) 

IMARK  =  I 
END  IF 

1030  CONTINUE 

C  IF  GREATER  THAN  LAST  POINT,  USE  LAST  POINT 
C  CENTERED  DIFFERENCE 
1040  CONTINUE 

IF  GMARK  £Q.  NEL)  THEN 
AMT  =  ELEC  (NEL) 

AMTDER  =  (FUEC  (IMARK)  -  ELEC  (IMARK  -  2))/(2.*(ELTIM  (IMARK) 
1  -  ELTIM  GMARK  -  1))) 

ELSE 

AMTDER  =  (ELEC  (IMARK  +  1)  -  ELEC  (IMARK  - 

1  1))/(2.*(ELTIM  GMARK) 

2  -  ELTIM  (IMARK  -  1))) 

END  IF 

C 

GO  TO  1020 
C 

C  IOPT=2  DERIVATIVE  OF  WF 
1010  CONTINUE 
C 

DO  1050  I  =  2.  NWF 
IF  (TIME  .EQ.  WFTIM  (I))  THEN 
AMT  =  WFCUM  G) 

IMARK  =  I 
GO  TO  1060 
END  IF 

IF  (TIME  .GT.  WFTIM  (I  -  1))  THEN 
AMT  =  WFCUM  G  -  1)  +  (WFCUM  (I)  -  WFCUM  (I  -  1))* 
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((TIME  -  WFTIM  a  -  l))/(WFnM  G)  -  WFTIM  (I  -  1))) 

IMARK  =  I 
END  IF 
1050  CONTINUE 

C  IF  GREATER  THAN  LAST  POINT,  USE  LAST  POINT 
C  CENTERED  DIFFERENCE 
1060  CONTINUE 

IF  (IMARK  JEQ.  NWF)  THEN 
AMT  =  WFCUM  (NWF) 

AMTDER  =  (WFCUM  (IMARK)  -  WFCUM  (IMARK  -  2))/(2.*(WFTIM  (IMARK) 
1  -  WFTIM  (IMARK  -  1))) 

ELSE 

AMTDER  =  (WFCUM  GMARK  +  1)  -  WFCUM  (IMARK  - 

1  1))/(2.*(WFTIM  (IMARK) 

2  -  WFTIM  (IMARK  -  1))) 

END  IF 

C 

1020  CONTINUE 
RETURN 
END 

SUBROUTINE  OUTPUT  (lOPT) 

C 

IMPLICIT  REAL  (A  -  H.  O  -  Z) 

PARAMETER  (NPT  =  150,  NPTT  =  50) 

COMMON  /START/  CHVOL,  AREA.  WFMASI,  RHOI, 

1  PROJI,  PROJW,  TOL.  TIMINC, 

2  TIMEND,  WAIR,  BM,  DERBM 

COMMON  /CURRENT/  TIME,  PRESCH,  TRAVPR, 

1  ELECEN,  TEMP.  RIMPET.  GAMMA, 

2  CO  VOL.  WFl,  WF2.  ELDEN,  SPMPR,  VELP,  BETA,  BASEPR,  CM, 

3  Z.  RHO,  ENCOMP 

COMMON  /OPTIONS/  IBETA,  ICOM 
COMMON  /DATA/  PRTIM  (1:NPT),  PRES  (1:NPT). 

1  TRTIM  (1  ;NPT).  TRAV  ( 1  :NPT). 

2  VEL  (1  :NPT).  ELTIM  (1  :NPT).  ELEC  (1  :NPT) 

COMMON  /OUT/  WFTIM  (NPT),  WFCUM  (NPT) 

COMMON  /NUMBER/  NPR,  NTR,  NEL,  NTH.  ITRVEL.  NWF 
COMMON  /NAME/  PRESIN,  TRAVIN,  ELECIN,  THERIN,  INITIN 
CHARACTER*20  PRESIN,  TRAVIN.  ELEQN,  THERIN,  INITIN 

C 

C  FILE  28  IS  OUTPUT  FOR  DETAILS  OF  WORKING  FLUID 
C  FILE  27  IS  ENERGY  BALANCE 
C  FILE  26  IS  ERROR  FILE  AND  MESSAGES 
IF  aOPT  .EQ.  1)  THEN 
WRITE  (*,  6000) 

OPEN  (26.  FILE  =  ’MESSAGES.OUT’) 

OPEN  (27.  FILE  =  ’ENERGY.OUT’) 

OPEN  (28.  FILE  =  ’WORKFL.OUT’) 


60 


OPEN  (29.  FILE  =  ’COMPRESS.OUT’) 

GO  TO  1000 
END  IF 

IF  aOPT  .EQ.  2)  GO  TO  1010 
IF  aOPT  .EQ.  3)  GO  TO  1020 
IF  aOPT  .EQ.  4)  GO  TO  1030 
C 

1010  CONTINUE 
WRITE  (*,  6010) 

WRITE  (*.  6020) 

WRITE  (27.  6030) 

WRITE  (28.  6040) 

WRITE  (29.  6050) 

GO  TO  1000 
C 

1020  CONTINIJE 

C  CHANGE  TO  MPa  AND  TO  JOULES 
TIMMS  =  TIME*1.0E+3 
PRESCH  =  PRESCH*1.0E-7 
BASEPR  =  BASEPR*1.0E-7 
SPMPR  =  SPMPR*1.0E-7 

WRITE  (*.  6060)  TIMMS.  PRESCH.  BASEPR.  TRAVPR.  VELP,  ELDEN,  WF2 

write  (26.  6070)  timms.  presch.  basepr 

rimpet  =  rimpet*1.0e-7 

enused  =  (riinpet/(gamma  -  l.))*wf2 

enkp=  .5*projw*velp*velp*1.0e-7 

alpha  =  wf2Av£cnasi 

cm  =  (wfmasi*(alpha  +  beta*(l.  -  alpha))  +  wair)/projw 

enkg  =  (l./6.)*cm*projw*velp*velp*1.0e-7 

vol  =  chvol  -  covol*wG  •  (wfmasi  -  wf2)/rho  +  travpr*area 

enigas  =  spmpr*vol/(gamma  -  1.) 

encom  =  (1.  -  z)* wfmasi  *encomp*1.0e -7 

if  (enused  .ne.  0.0) 

1  enerr  =  ((enused  -  enkp  -  enkg  -  enigas  -  cncom)/enused)*l(X). 
write  (27.  6080)  timms.  enused.  enkp,  enkg,  enigas,  encom,  enerr 
write  (29,  6090)  timms,  rho,  encomp,  encom 
C 

GO  TO  1000 
C 

C  BF  IOPT=l,  FIND  DER  OF  ENERGY 
C  IF  IOPT=2,  FIND  DER  OF  WF 
1030  CONTINUE 
TIME  =  0.0 

TIMEND  =  WFTIM  (NWF) 

TIM2  =  0.0 
EL  =  0.0 
WF  =  0.0 
ELDEN  =  0.0 
ELDER  =  0.0 
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WFDER  =  0.0 

WRITE  (28,  6100)  TIM2,  EL.  WF,  ELDEN,  ELDER.  WFDER 
DO  1040  I  =  1.  MPT 
TIME  =  TIME  +  TIMING 
IF  (TIME  1£.  TIMEND)  THEN 
CALL  DER  (TIME,  1,  EL,  ELDER) 

CALL  DER  (TIME,  2.  WF,  WFDER) 

IF  (WF  JSTE.  0.0)  ELDEN  =  EL/WF 
TIM2  =  TIME*l.E+3 
C  ELDEN  IN  MJ/s,  WFDER  IN  kg/s 
ELDER  =  ELDER*.001 
WFDER  =  WFDER*.001 

WRITE  (28.  6100)  TIM2,  EL.  WF.  ELDEN.  ELDER.  WFDER 
ELSE 

GO  TO  1000 
END  IF 

1040  CONTINUE 
C 

6000  FORMAT(/.*  Inverse  Electrothennal  Gun  CodeV/) 

6010  FORMAT(lX.’  TIME  *.2X,’CHAM  PR*.2X. 

1  ’  BAS*.  PR  ’.aX.’PROJ  TRAV*,2X.’  PROJ  VEL’.5X.’EL  EN  DEN’, 

2  2X.’CUM  WF*) 

6020  FORMAT(2X.’  (ms)  ’.2x.2x.’(MPa)’.3x.2x.lx,*(MPa)’,lx.2x. 

1  2x.’(cm)*.3x.2x.3x.’(cm/s)’.2x.2x.2x.’(MJ/kg)*,lX, 

2  2X.lX.*(g)*y) 

6060FORMAT(1XJ='6.2.2X.F7.2,2X,F10.5,2XJ9.4.2X,F10.2,2X4^.3.2XJ^10.5) 

6070  format(lx.f6.2.2x,fl5.7.2x,fl5.7) 

6030  fonnat(lx.*ENERGY  BALANCE’^/, 

1  lx,’time’,7x,’enused’,10x,’KE  proj’,9x,’KE  gas’,8x, 

2  *Intgas’.5x.’En  comp’,5x,’Error 

3  /.lx.*(ms)’.7x.’ (J) ’.llx.*  (J)  ’.9x.’(J)  ’.lOx,’  (J)*. 

4  lOx,’  (J)*) 

6080  fonnat(lx.fl0.3.1x.fl5.2.1x,fl0.2.1x.fl0.2.1x.fl2.2.1XJ^10.2.1x, 

1  F6.1) 

6040  FORMAT(lx.’WORiaNG  FLUID  AND  ELECTRICAL  ENERGY  SUMMARY’7/, 

1  lx.lx.”nME’,lx,3x.’EL  EN’ Jx.5x.’WF  CUMULATIVE’.lX.lx, 

2  ’EL  EN  DENSITY’,lx.’ELEC  DER’,lX,2x,’WF  DER’, 

3  /,lx,lX.’(ms)’.lX.3x.’(kJ)’.3X,l  lX.’(G)’.6x.4X.’(MJ/kg)’.5X. 

4  2x,’(MJ/s)’,2X,3x.’(kg/s)’) 

6100FORMAT(lx.f6.3.1x,fl0.2.1x.fl5.5.1x.n5.5,lx.n0.2,lx.n0.2) 

6050  FORMAT(lx,’timms(ms),rho(g/cc),encomp(J/g),cncom(J)’) 

6090  FORMATd x,f6.3.1  x,f8.5,l x,f20.5,l x.flO.5) 

C 

1000  CONTINUE 
RETURN 
END 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  uiKlertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-TR-3391 _ Date  of  Report  September  1992 _ 

2.  Date  Report  Received _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ _ _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so.  please  elaborate. _ 


.  6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct  Address 
in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


Gty,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 
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OFRCIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

BUST  CUBS  Md  01101,  AP6.  MD 

Postage  mHI  be  paid  by  addressee. 
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